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1. Introduction 

It is well known (see, for example, [|T], |2|, |3|) that in a nonlinear medium with the Kerr 
nonlinearity, i.e. with cubic nonlinearity, self-phase modulation (SPM) phenomenon 
occurs which gives rise to quadrature squeezing of fluctuations with conservation of 
the photon statistics. This effect is observed when the one-mode radiation propagates 
in the Kerr medium. If we deal with the propagation of the two-mode radiation in 
such a medium, the situation becomes more complicated. Besides the SPM effect, the 
phenomena of cross-phase modulation (XPM) and parametric interactions can, generally 
speaking, take place. However, the efficiency of parametric frequency-multiplication 
processes depends on the phase mismatch of the interacting waves. The parametric 
frequency conversion can be neglected in the case of the strong phase mismatch. It is 
worthy noting that the influence of the parametric energy exchange at the two-mode 
interaction on the generation of polarization-squeezed light in the Kerr medium has 
been considered in |4[]. Both SPM and XPM effects do not depend on the phase 
mismatch, consequently, they can play a dominant role when one investigates the 
nonlinear propagation of two modes with orthogonal polarizations and/or different 
frequencies in the Kerr medium. To the best of our knowledge, such analysis has been 
carried out only for the case of monochromatic modes (see review ||). 

In order to develop a correct quantum theory of the nonlinear propagation of 
the pulse radiation, it is necessary to take into account a finite response time of the 
nonlinearity. There are two consistent approaches to solving this problem. In |], 0] 
the Kerr nonlinearity has been treated as a Raman -like nonlinearity. The authors 
H |7j have taken into account quantum and thermal noises as a fluctuating addition 
to the relaxation nonlinearity in the interaction Hamiltonian. In the other approach 
developed in || |9|, [10[| , the interaction Hamiltonian || or the interaction momentum 
operator ||, [lOf] , which contains only a response function of the nonlinearity, have been 
used. In such quantum description of the SPM of ultra-short pulses (USPs), one cannot 
add thermal noise terms to the interaction Hamiltonian and momentum operator in 
order to satisfy the canonical commutation relation for the annihilation and creation 
Bose-operators. 

In this paper, we present the quantum theory based on the interaction momentum 



operator [[J for the case of nonlinear propagation of two ultrashort light pulses [10 
In reality, here we will mainly analyze two combined (SPM and XPM) effects. The 
first one is responsible for the generation of the squeezed state of light and the second 
one (as we demonstrate below) can be employed to control the squeezing process. This 
combined phenomenon can be shortly called "the SPM-XPM of USPs." Since we deal 
with the analysis of the combined phenomena, algebra of time dependent Bose-operators 
developed earlier in || [IT] has to be subsequently extended. 
The paper is organized as follows. 

In section |2] the quantum equations for the SPM-XPM of USPs are introduced 
when the nonlinearity of medium has relaxation behaviour of an electronic Kerr type. 
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Algebra of time- dependent Bose-operators is extended in section |3| and used to estimate 
statistical parameters of the USP quadrature components. In section [| the correlation 
functions of the quadrature components are obtained. Then in section |5] the spectrum 
of quadrature-squeezed component of the USP under investigation is studied. Section ^| 
closes the paper with some concluding remarks. 



2. Quantum equations for simultaneous SPM and XPM of light pulses 

The conventional way to derive the quantum equation of SPM is based on the interaction 
Hamiltonian, with the quantum equation describing the time evolution of annihilation 
or creation Bose-operators being derived. A transition to the space evolution is usually 
provided by the replacement t — > z/u, where z is the distance passed and u is the pulse's 
group velocity in a nonlinear medium. This approach seems to be fairly appropriate for 
the case of single-mode radiation. Generally speaking, both t and z are present in the 
analytical description of nonlinear propagation of the pulse. Therefore, the momentum 
operator connected with the spatial evolution of the pulse field should be used instead 
of the interaction Hamiltonian [12| in order to obtain an equation for the annihilation 
Bose-operator. This approach for the case of pulse SPM in nonlinear medium with 
inertial electronic Kerr nonlinearity has been successfully used in ||. 

We start with the analysis of SPM-XPM of two coherent USPs in a noninertial 
nonlinear medium. 



2.1. Quantum SPM-XPM equation in noninertial nonlinear medium 

In a noninertial nonlinear medium, the combined SPM-XPM effect of USPs is described 
by making use of the momentum operator G(z) JK| 

G(z)=j2Ggl(z) + G xpm (z), (1) 
i=i 

where index j denotes the pulse's number (j = 1,2), and the momentum operators of 
SPM and XPM effects are 

Ggl(z)=hpJ N[nJM]cft, (2) 

J — oo 
~ roc 

G xpm (z) = hf3 hi(t, z)n 2 (t, z)dt. (3) 

J — oo 

Here h is Planck's constant, N is the operator of normal ordering, the factors /3j and j3 are 
defined by the electronic Kerr nonlinearity n 2 of the medium (for example, (3j = C 2 7j/2; 
C = (hu /2V) 1 ^ 2 , 7* = k n 2 (j)/2n |J), and fij(t,z) = Aj~(t, z)Aj(t, z) is the photon 
number "density" operator for the j-th pulse in a given cross-section z at the time 
moment t. The operator Aj(t,z) [Aj(t,z)] is the photon creation (annihilation) Bose 
operator. It is the slowly varying operator of the negative (positive) frequency part of 
the electric field strength operator The XPM momentum operator defined by (|||) 
takes into account the action of the first pulse on the second one and vice versa. 
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iPih^t, z) + i(3h 2 (t, z) A 1 (t,z)=0. (5) 



In the Heisenberg representation, the space evolution for the time-dependent Bose- 
operator Aj(t, z) for the j-th pulse is given by the equation (see fl~2|j ) 

-ih^^=[A j (t,z),G(z)]. (4) 

In accordance with (|l|), the space evolution equation, for example, of the operator 
Ai(t, z) reads 

dA 1 {t,z) 
dz 

The spatial evolution equation of the operator A2(t,z) can be easily obtained by 
permuting the indices 1 •*->• 2. Thus, we get the system of four coupled equations keeping 
in mind two equations for the creation operator. One points out that Eq. (|5p is written 
in the moving coordinate system: z = z and t = t —z/u, where t' is the running time, u 
is the pulse's group velocity in the nonlinear medium. We assume that one can neglect 
the difference in group velocities of the pulses and do not consider the pulse spreading 
due to the medium dispersion, that corresponds to first approximation of the dispersion 
theory. 

In view of Eq. fl5|), one can show that the operator hj(t, z) does not change in the 
nonlinear medium: 

hj(t, z) = hj(t, z = 0) = n 0J (t), (6) 

where z = corresponds to the input of the nonlinear medium. Equation (||) means 
that the photon statistics does not change in the medium. 

Solving Eq. (||) and its Hermitian conjugate, one has for the annihilation and 
creation photon Bose-operators 

ii(t,*) =e^ 1 ^.»( t )i 0i i(i) l (7) 
At(t, z) = A+ tl (t) e -nmo,iW-nno, 2 (t) ) ( 8 ) 

where 71 = f3\z and 7 = j3z. By permuting the indices 1 <-» 2, one gets the Bose- 
operators for another light pulse. 

The preservation of the canonical structure of quantum theory requires that the 
annihilation and creation Bose operators Aj(t, z) and Aj(t, z) satisfy the commutation 
relations 

'Ajfa, z), A+(t 2 , z)] = 8 jk 6(t 2 - h) (9) 

for an arbitrary distance z in the medium, where 8^ is the Kroneker delta-symbol, 
j,k = l,2. 

There are some peculiarities in the quantum description of the combined SPM- 
XPM effect in a noninertial nonlinear medium. Firstly, the solutions (^), (|8|) do not 
permit to verify the commutation relation (|9|). Secondly, the reduction of the expressions 
e mno,j(t)^ e *7flo,j-(t) to the normally-ordered form is accompanied by the appearance of 
a nonintegrable singularity (see || [TT|). These circumstances do not appear in the 
quantum theory of combined SPM-XPM effect which takes into account the relaxation 
behaviour of the nonlinearity. 
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2.2. Quantum equation of pulse SPM-XPM in inertia! nonlinear medium 

The momentum operator for SPM of USP in a medium with the electronic Kerr 
nonlinearity has been introduced in ||. It incorporates the function of a nonlinear 
response in its structure. In the model considered, a contribution of the electronic 
Kerr nonlinearity decreases exponentially. Thus, we rewrite the momentum operator of 
SPM-XPM effect G(z) introduced according to (0) (see 0) as follows 

G®{z) = hPj r dt f H(t - ti)N[hj(t, z)]dt u (10) 

A ~ /'OO rt 

G xpm (z)=hf3 dt H(t-t 1 )[hi(t,z)h 2 (ti,z)+hi(t 1 ,z)n 2 (t,z)]dti, (11) 

J —oo J —oo 

where H(t) is the function of nonlinear response, asymmetrically defined in order to 
satisfy the condition imposed by the causality principle: H(t) ^ at t > and 
H(t) — at t < 0. The term in the second integral in ( |TUD should be interpreted 
as a generalized force acting in the cross-section z of the medium which at the time 
moment t depends only on the previous time moments. The similar term in the second 
integral in ([TT1) should be interpreted as a sum of two generalized forces acting between 
pulses in the same cross-section z of a medium which also depend on the previous time 
moments. Therefore, the causality principle in the Hermitian operators ([T0|) and (p|) 
is not violated. 

In the case of the nonlinearity of an electronic origin, the nonlinear response function 
can be introduced as follows: 

H{t) = (1/tv) exp (-t/r r ), (t > 0) (12) 

and H(t) = at t < 0. Indeed, if a single USP propagates through the Kerr medium, 
then the evolution of nonlinear addition An n i to the refractive index, which is associated 
with the SPM effect, is given by the equation (see, for example, [§, [U], |T3|1 ) 

dAn nl {t,z) 1 

+ An n i{t,z) = -n 2 (j)n j (t,z), (13) 



dt ' — 2 
which has the solution 

1 /•* 

An nl (t, z) = -n 2 (j) / H(t -ti)n j (t 1 ,z)dt 1 . (14) 

Note that the nonlinear response function (^) appears in the absence of one- and 
two-photon and Raman resonances |13|]. Therefore, the approach developed is valid 
when each pulse frequency is off-resonance and the pulse duration r p is much larger 
than the relaxation time r r . If the USP propagates in a fused silica- fiber, then about 
80% of the Kerr effect is due to the electronic motion occurring on ~ 1 fs time scales 



and only 20% of the Kerr effect is attributable to the Raman oscillators | i~4| . Thus, 
our model corresponds to the case of the Kerr effect mainly produced by the electronic 
motion. 

Since in quantum theory the relaxation behaviour is connected with the so-called 



thermal reservoir, the operator evolution equation (13) must, in general case, contain 
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a source of thermal noise besides a relaxation term. Then we can write the following 
expression for the nonlinear addition |6|. [/J 

1 ft 

Ah nl (t, z) = -n 2 (j) / H(t -t 1 )n j (ti,z)dt 1 + m(t,z), (15) 



2 

where the Hermitian operator rh(t, z) takes into account thermal fluctuations of 
An n i(t,z) in the absence of light field. For the nonlinearity of electronic origin 
an expression likes fllBD can be obtained from Duffing-type equation fl3| . Here it 
is important to note that the average value of m(t, z) operator is equal to zero; 
(m(t,z)) = (see also 0). Hence, one can consider Eq. fll~4]) as the result of averaging 
Eq. (|15D over the thermal fluctuations. In connection with that, expressions fllCf ), ([ll]) 
can be truly considered as ones averaged over the thermal fluctuations of nonlinearity. 

In the present paper we neglect the thermal fluctuations since we are interested 
in the nonlinear phase fluctuations caused by the quantum USP ones. Nevertheless, 
the commutation relations @ are proved to be valued in this simplified approach (see 
below). However, it is not difficult to generalize the theory of the SPM and XPM 
phenomena allowing for thermal noise in the Kerr medium with electronic nonlinearity. 

Using (J^j it is easy to prove that the operator hj (t, z) commutes with the the 
momentum operator G(z), 

[6(z),hj(t,z)] = [G{z),fioAt)]=0, (16) 

that is, the photon number operator fij(t, z) remains unchanged in the nonlinear medium 
[see also Eq. ([])]. 

Taking into account fllOD and (|TT|), one obtains from (f|) the quantum equation for 
combined SPM-XPM effect, for instance, for the pulse with index 1 

dA^z) 
dz 

where 



{i/Mno,i(*)] +i/fy[no,2(*)]} = 0, (17) 



#o,j(*)] = / H(h) [n 0J (t - h) + n 0ij (t + ti)] dt x . (18) 

The appearance of the second term in ( |18D related to n 0t j(t + ti) in the quantum 
description has been already discussed in where it has been assumed that this term 
can be connected with the vacuum fluctuations which are present even in the absence 
of a pulse. It should be also pointed out that expression ( |T7D is an intermediate result 
written in the moving coordinate system. 

Solving the space evolution equation QT7| ) for the annihilation Bose-operator and 
its Hermitian conjugate we obtain 

A x (t, z) = e Hi9N,i(t)]+m[no, 2 (t)] A 0tl (t), (19) 
A+(t, z) = A+ tl (t) e -higlno,i(t)}-iyg[no,2(t)} ^ (20) 

One can obtain similar expressions for the Bose-operators of the second pulse by 
permuting indices 1 •*->• 2. It is convenient to rewrite expression (|IBD as follows 

roo 

q[n 0J (t)} = / MtiK^t - t x )dt x , [h(t) = H(\t\)}. (21) 
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If we consider n j to be time independent, then ( Jl9| ) and ( p0[ ) describe the case of 
monochromatic modes (see, for example, |5|, |15[). If in ([Tj^) and (^) the response 
function Hit) = S(t), then the description of the noninertial nonlinear media (0) and (|8|) 
can be produced. Note that in the quantum description the structure of the nonlinear 
response ( |21~D is similar to the one of a linear response of a medium in second-order 
approximation of the dispersion theory (see || [T3| ) . 

To estimate the statistical characteristics of pulses at the output of the nonlinear 
medium, we need to calculate the average values of the operators' moments. As it is 
well known, they can be found if the operator expressions are given in the normally 
ordered form. The use of solutions (pl|) and (|20|) involves the development of a special 
mathematical technique. Bellow some elements of algebra of time-dependent Bose- 
operators developed in |§ are extended. 



3. Algebra of time-dependent Bose-operators 

For the quantum analysis of the pulse's SPM the correspondent algebra of time- 
dependent Bose-operators has been developed in |J. In the case where the XPM effect 
is present in addition to the SPM effect, we need to improve the mentioned algebra. It 
is convenient to introduce the operators 

d,(t) = z 7j q[fioj(jt)], Oft) = Hq[n 0>j (t)l (22) 
and their Hermitian conjugates 

6 jit) = -0,(t), O^t) = -d,(t), (23) 

in order to simplify further calculations. Thus, Eqs. ([19]), (20) may be represented in 
the operator form: 

A l (t,z) =e 6 ^ +d ^A , 1 (t), (24) 

At(t,z) = A+^e 6 ?®* 6 *®, (25) 

where the operator 0\ it) is responsible for the SPM of the j-th pulse and the operator 

Oii$) corresponds to the XPM effect due to the action of the second (control) pulse on 
the first pulse (under investigation). 

We suppose that the initial pulses are in coherent states, the j-th pulse's operator 
A)j(t) acting only on the state vector \a j(t)) within the associated Hilbert space Tif 

A ,j(t)\a ,j(t)) = a ,j(t)\a ,j(t)), 

where ctoj(t) is the eigenvalue of the operator. The factorization of the states of two 
different sub-Hilbert spaces Hi and 7i 2 takes place. The summarized quantum state of 
two pulses is described by the vector 

\a (t)) = |oo,i(t)) ® |ao,2(t)) 
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in the global Hilbert space Ti = Hi <S> 7i 2 - Since the coherent quantum states of each 
pulse occupy a distinct sub-Hilbert space for any two arbitrary chosen time moments ti 
and t 2 one can write : 

[hi{t u z),h 2 {t 2 ,z)\ = [fto,i(*i),no,2(*2)] =0 ' ( 26 ) 
from which we obtain 

[A 0>1 (ti), e 6 ^}= 0, [Aoflih), e 6 ^}= 0, (27) 
e 6 ^]= 0, [i+^tO, e^] = 0. (28) 

Calculating ( pT| ) and ( pB"D it is easy to prove that the following relationships take place: 

e 6 2 ( t2 )] = [e 6+ (tl)) e o+(fe) ]= 0> (29) 

[ e o l(tl)) e o 2 ( t2 )j = [e dat l); e o 2 (fe)j = (30) 

It should be pointed out that in the next sections the averaging operations are made 
over the total quantum state \ao(t)). 

3.1. Operator permutation relations 

For each pulse, the following operator permutation relations hold (see also |§): 

Aoj^Ojih) = [6 3 {t 2 ) + Vj(t 2 - tOAjfo), (31) 
6^^)^+ (t a )=il+(t a )[6 i (ti)+2? i (t3-ti)], (32) 

where 

Vjih-h) =iljH t 2-ti)- 

Besides, T>j(t 2 — ti) = ^(ti — ^2) since /i(t) is an even function of t [see In view 

of the mathematical induction, it is possible to demonstrate the validity of the formulae 
(m G N): 

Aoj^Of (t 2 ) = [dj(t 2 ) + Vj(t 2 - hTA^ih), (33) 
Of (^A^ih) = A+^idjih) + Vj(t 2 - h)] m . (34) 

Expanding e°^ and e°i ® in Taylor series, we can get the operator permutation 
relations which play an important role while estimating the statistical characteristics 
of the pulse investigated. By using (|33|) and (|34]) we finally arrive at 

i 0j (ti)e^ ( * 2) = e^+^-^ioj^i), (35) 
e^A+jih) = A+^e 6 ^^^. (36) 

Other permutation relations can be obtained by Hermitian conjugation. Making use of 
the relations (|35|), (|36|) and Eqs. (|29"D, fl3~0|) one can show that the canonical commutation 
relation (|9]) for the operators Aj(t, z) and Af(t, z) is exactly fulfilled. Proceeding in the 
same manner we can prove the validity of the following relations 

[A j {t 1 ,z),A j {t 2 ,z)}=0, (37) 
[hj(ti,z), hj(t2,z)] =0. (38) 
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Let j=l. Then 

[ii(ti, z), ii(t 2 , z)\ = Axih, z)A 1 (t 2 , z) - A^h, z)A 1 (t 1 , z) 

= e o 1 (* 1 )+<3 2 (t 1 )i 0]1 ( tl ) e o 1 (t 2 )+(3 2 ( t2 )^ i(t2) 

_ e 6i(fo)+d 2 (fa)^ i(t2)e 6 a ( tl ) + 6 2 ( tl )^ 0)i(fi) 

= e 6 1 (tx)+6 1 ( t2 )+d 2 ( tl ) + d 2 (t 2 ) e ^(* 2 - tl )^ oi(ti) ^ 0)i ^ 2) 

_ e 6 1 ( tl )+6 1 ( i2 )+(3 2 ( tl ) + (3 2 ( t2 ) e © 1 ( tl - t2 )^ i(ti)y | 0ji(t2) 

= 0. (39) 
We used above the permutation relation (^). Let us verify our initial statement (d). 
hx(t,z) = A+(t,z)Ax(t,z) 

= i+ 1 (t)e 6 fW + ^(*) e ^(*)+^Wi 0il (t) 

= i+ 1 (t)e 6l(i)+ ^ w+d2W+ ^ (t) i 0il (t) 
= i+ 1 (t)i ,iW 

= n ,i(t). (40) 

Hence we have 

[ni(*i,2:),ni(t2,£)] = K,i(^i), ^0,1(^2)] 

= ^0,1^1)^0,1^2) - ^0,1(^2)^0,1^1) 
= i+ 1 (i 1 )i 0il (t 1 )i+ 1 (t 2 )io,i(t 2 ) 

-i^^)^,!^)^!^)^,!^) 

= i+^tOi+^ta)^!^)^!^) + S(t 2 - t^A+^Ao^h) 

= 0. (41) 

In (f4T|) we used the fact that for any functions f(x), g(x) we have S(x 2 — Xi) f (xi)g(x 2 ) = 
S(xi — x 2 )f(x 2 )g(xi), where x\ and x 2 are two arbitrary values in which the functions 
are defined. 

3.2. Normal ordering 

As it was stated above, another important issue is presenting the operators Aj(t, z) and 
A~j~(t,z) in the normally ordered form. The normal ordering formulated in the time- 
representation allows one to estimate the means of the Bose-operators over the initial 
coherent states. As a consequence, the operator e° j ^ takes the form (see P, |§): 

eAW = Nexp ||°° [e H{e) - l] n 0>j (t - 9r r )de\ , (42) 

where 9 = t/r r , 7i(6) = i^h{9) and h(6) = r r h(6T r ). The operators in the integral 
in (j|^) should be understood as the c- numbers. Thus, averaging the e° J W over the 
coherent state results 



e 



Oj(*)\ 



00 



) = exp I J [e nw - lj n ,j(t - 9r r )d6j (43) 
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where 

no A*) = ("oj(*)) = l«oj(t)| 2 
is the average photon number density of the j-th pulse at the input of nonlinear medium 
z — 0. The factorization of the initial coherent states allows one to write 

(a (t)\e 6 ^ +d ^\a (t)} = (e^)(e da % (44) 

where 

(e dl W> = <a , a (t)|e^W|a b , 1 (t)>, (45) 

(e d * (t) ) = (a , 2 (t)\e d ^\a , 2 (t)). (46) 

3.3. Means of Bose-operators and their combinations 

In majority of experimental cases, the parameter 7 C 1 and due to this we can 
decompose the term in integral in ( f43"D and truncate the decomposition by terms of 
the order of 7 2 (see 0). As a result, we get 

/ e 6i(t)\ e ^i(*)-M*) ; / e 6+(t)\ e -^i(t)-M*) ; (47) 

/g0 2 (t)^ _ e # 2 (*)-A«2(t)^ ( e 2 (*)^ _ e -i^ 2 (t)-A2(t) ; (43) 



where 



0i(£) = iUo.i / M% 2 (i - 0r r )tZ0, (49) 



2 



-00 
00 



2 (t) = ^0o,2 / h{6)r 2 2 {t - 9r r )d9, (50) 
Z J —00 

Pi(*) = oPo,i / h 2 (9)rl(t-9T r )de, (51) 



2' 



—00 
00 



= ^0,2 / ^(^)r 2 2 (t - 9r r )d6. (52) 
Here we also introduce the parameters 

0o,i = 27^0,1(0), fiQ tl = 7 2 "-o,i(°) = 7i0o,i/2, (53) 

00,2 = 27^0,2(0), fi ,2 = 7V),2(0) = 700,2/2, (54) 

with rj(t) being the envelope of j-th pulse so that oco,j(t) = ao,j(0)rj(t) and rj(0) = 1. 
Let us denote for simplicity n Q j = n 0j (0). The parameters 0i(t), f^x(t) are certainly 
connected with SPM of the investigated pulse and the parameters 02 (t), fait) are 
connected with the XPM of pulses. The parameters 4>\{t) and 02 if) have the physical 
meaning of the nonlinear phase additions caused by the SPM and XPM, respectively. 
Here 271 represents the nonlinear phase shift per one photon for the investigated pulse. 

A particular interest is connected with the estimation of the average values of 
combinations of the exponential Bose-operators over the coherent state. Using the 
procedure of normal ordering we have the following formulae: 

/ e 6x(tl)+(5l(t2)\ = e i[M*i)+^i(*2)]-Mi(*i>*2)-£i(ii,t2) (55) 

( e dt(h)+di(t2)} = e il-Mti)+Mt2)]-Mh,t2)+iCi(ti,t2)^ ( 56 ) 
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3 O a (ti)+0 2 (t 2 ) 



} 



J [-02 (t 1 )+<fe (*2 )] -M2 (*1 ,<2)+K 2 it i ,t 2 ) 
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(57) 
(58) 



where 



l^i(ti,t 2 ) = /ii(ti) + /ui(t 2 ), At2(*i, *2) = £2(^1) + 
and /Ci(ti,t2), ^2(^1,^2) are the time correlators: 

/Ci(*i, £ 2 ) = /i ,i / - 0T r )/i(t 2 - 9T r )rl(6T r )d9, 



JC 2 (t 1 ,t 2 



-00 
00 



- 9r r )h(t 2 - 9r r )rl(9r r )d9. 



(59) 
(60) 



As it was mentioned above, our theory makes use of the assumption that t p ^> r r . 
Therefore one can simplify expressions (^9|) — ( 



|) and (]59f) , ( pQ) eliminating 7j(£ — 6>r r ) 
and rj{9r r ) from the integrand in the particular points: 6*r r = in (f49j)-(|52l), and 
9r r = t\ + t/2 in fl59| ), (|60|), where r = t 2 — t 1 . Taking into account the response 
function for the electronic Kerr nonlinearity fll2|) for the expressions mentioned, we 
obtain 



Mi) =0o,ir?(t), 
=/i ,ir?(*)/2 

and the correlators 

/Ci(ti,t 2 ) =/i ,ir 1 2 (ti + r/2)( ? (r), 
M*i,* 2 ) =/io,2r 2 2 (ti + r/2)( ? (r). 
The function g(r) has the form 

1 

T, 



02 (t) = 0o,2 r£(t), 

At2(*) = /io,2^2^)/ 2 > 



*(t) 

















TV , 





(61) 
(62) 

(63) 
(64) 

(65) 



The results obtained permit to investigate the statistical properties of the USP in a 
medium with the Kerr nonlinearity. 



4. Correlation functions of quadrature components 

Here we restrict our analysis by studying the quadrature components which are defined 
by the expressions: 

1 r 



X,(t, z) 
Yi(t,z) 



2i 



A 1 (t,z)+A+(t,z) 

A 1 (t,z)-A+(t,z) 



(66) 
(67) 



As it was mentioned above, in the nonlinear media under consideration the photon 
statistics of each pulse remains unchanged. 
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The mean values of the operators Ai(t, z) and Af(t, z) for the case of initial pulses 
being in the coherent state are 

(A 1 (t,z)) =a , 1 (t)(e 6 ^)(e d ^), (68) 

(At(t,z))=a* Qtl (t)(e d ^)(e^^). (69) 

Taking into account Eqs. (f4T|), ( |48| ) and the fact that aoj(t) = \aoj{t)\e tlfj ^\ for mean 
values of quadratures we get 

(Mt>z)) = |a 0l i(*)|e^ l(t) - A2(t) cos [$ x (t) + &(*)], (70) 
= laoiOle-^^^^ sin + &(*)], (71) 

where = 0i(t) + and <^i(i) is the linear phase of the investigated pulse. The 

exponential terms in (f70|) and (]TTf) are caused by quantum effects at the SPM and XPM 
(there are no such terms in the classical theory). From (f70| ) and (]71|) one can conclude 
that the changes of quadratures in time are quasi-statically connected with the changes 
in envelopes of both pulses [<f>i(t) = <f) 0) irf(t) = 271^0,1 (t), 4>2(t) = 4>o,2 r 2 (t) = 27^0,2 (£)]■ 
Remind that t is time in the coordinate system moving with the group velocity. In other 
words, the causality principle is satisfied for the parameters observed. 
We define now the correlation functions of quadrature components 

R Xl {t, t + r) = (Xx(t, z)Xi(t + t,z))- (!i(t, z)){X x {t + r,z)), (72) 
R Yl (t, t + r) = (Yx(t, z)Yi(t + r, z )> - (Yi(t, z)) (Y\(t + r, z)>. (73) 

The analysis of the correlation functions requires the evaluation of the correlators 
C Xl (t u t 2 ) = (X^Mh)) and C Yl (t u t 2 ) = (^(tOYiCta)). 

Using the permutation relations fl3~5f ), (|36|) and Eqs. (|55| ) - (|58l) we obtain 

Cx^ti, h) = - A 5{t 2 ~ti) + \ |oo,i(*i)| Kite)! e-w.^' 1 ^) 

X {e~ Al ' 2(il ' i2) cos [$ 1)2 (t0 + $ li2 (t 2 ) + 7l fc(t 2 - h)\ 
+ e A 1 , 2 (ti,* a ) C0S [$ li2 (t 1 ) -<£ 1>2 (t 2 )]}, (74) 
1 

4 ( 



t 2 ) = 7<5(t 2 - *i) - J |«o,i(ti)| |«o,ife)|e-^ 2(tl ' i2) 



x < e 



e 



Ai, 2 (ti,t 2 ) 



2) cos[$i,2(*i)-$i,2(i 2 )]}, (75) 



where we introduced for simplicity the notations 

Ai,2 (ti, £2) = HiihM) + fcihM), 
Ai, 2 (*i, * 2 ) = Ai(*i, * 2 ) + Mh,t 2 ), 

Ai(ti,t 2 ) =Aii(ti,t 2 )^(t2-ti), (76) 

A 2 (ti,t 2 ) = /2 2 (<l,*2)^(*2 -*i)> 

$ 1)2 (t) = $i(t) 
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As a result, for the correlation functions of quadrature components we finally get: 

Rxi (t, t + r) = -{8(t) - <j)i{t)h{r) sin 2$ lj2 (t) 

+ [<f>l{t) + Mt)Mt)]g(r) sin 2 <§ 1>2 (t)}, (77) 

R Yl (t,t + r) = -{<f(r) + sin 2$ 1)2 (*) 

+ [0 2 (t) + 0i(t)02(t)]^(r) cos 2 $ li2 (t)}. (78) 
To obtain (ff7|) and (0), the 7 <C 1 and r r <C r p approximations have been used. 

5. Spectrum of quantum fluctuations of quadrature components 

The spectral densities of quantum fluctuations of the quadrature components are defined 
by the expression: 

/oo 
Rx^&t + ry^dT. (79) 
-00 

Taking into account a slow change of the envelope during the relaxation time, one arrives 
at: 

S Xl (u,t) = i{l-20 1 (t)L(w)sin2$ 1>2 (t) + 4[^ (80) 

S Yl (u,t) =^{l + 20 1 (t)LHsin2<l li2 (t)+4[0 2 (t)+0 1 (t)0 2 (t)]L 2 Hcos 2 <l li2 (t)}, (81) 

where L{uS) = [1 + (u;r r ) 2 ] _1 and 4>j(t) = 2jn Q j(t). 

From fl80|) and ( Ell ) it follows that in the case 2 (t) = 0, i.e. in the absence of XPM 
effect, the known result at the SPM effect for the investigated pulse can be obtained 
0, [|. At the SPM process, one can control the spectrum by choosing the phase of 
the initial investigated light pulse to be optimal for the chosen frequency 0, [TTfl . The 
presence of the XPM effect adds new terms in the multiplier and the phase in expressions 
(p0|) and (|8~ID in comparison with the ones for the SPM effect only. This circumstance 
gives us another possibility to control the fluctuation spectrum of the investigated pulse 
by varying these terms. 

From (|80D and (^l|) it also follows that the choice of the phase $i, 2 (t) determines 
the level of quantum fluctuations to be lower or higher than the shot-noise level 
Sx-l{w) = Syi(v) — 1/4, corresponding to the coherent state of the initial pulse. In 
agreement with the Heisenberg uncertainty relation, the spectrum of the Y^-quadrature 
is shifted by a phase it/ 2 in comparison with the spectrum of the Xx-quadrature. 

If one chooses the optimal phase of the initial investigated pulse 

= r rctan ( A ) " m ~ h[t) m 

for the reduced frequency f2 = ujqt Ti then spectral densities fl80D , fl8lD take the form 
S° Xl (no,t) = \{l - 2L(n M(t) + < 2 (t)L 2 (fi )] 1/2 +2^ j2 (t)L(fio)}, (83) 
S Yl (Q ,t) = i{l + 2L(fi o )[0l(t) + < 2 (t)^ 2 (fio)] 1/2 +2^ i2 (t)L(fi o )}, (84) 
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Figure 1. Fluctuation spectrum of the quadrature-squeezed component of the 
investigated pulse for the case of optimal initial phase <po,i(i) chosen for Slo = and 
for different intensities of the control pulse no, 2/^0,1 — ( a )> 2 (b), 3 (c), 5 (d), 8 (e). 
Curves are calculated for the time moment t — and 71 = 72 = 27, 0q,i = 271 no. 1 = 2. 

where <f)* 1)2 (t) = 4>\{t) + fait) fait). 
At any frequency Q, we have 

s Xl (n,t) = s° Xl (n ,t) + ±[L(n)-L(n )}{<pi 2 (t)[L(n) + L(n )} 

- + 4>Ut)L(n )[L(n) + L(n ))M(t) + < 2 (t)L 2 (Q Q )]- 1/2 }, (85) 
s Yl (n,t) = s° Yi (n ,t) + ^[L(n) - L{n )]{(f>i 2 (t)[L(n) + L{n )) 

+ {<Pl(t) + <Pi 2 (t)L(Q )[L(Q) + L(Q )]M(t) + < 2 (t)L 2 (fi )]- 1/2 }. (86) 

From (^) and flSB]) it follows that the change fait) = 27710,2 it) provides a possibility 
for controlling the spectrum of quadrature-component fluctuations. Certainly, this 
change can be realized by varying the control pulse intensity fiQ 2 it) at the optimal 
phase of the investigated pulse. 

The spectra of the investigated pulse at a fixed optimal initial phase 990,1 it) (chosen 
at Qq = 0, 0.5, 0.7 and the time moment t — 0) for different values of the photon 
numbers of the control pulse are displayed in Figs. ^ 0, and Q respectively. 

From Figs. it follows that the change of the intensity of one light pulse gives a 
method for the control of squeezed spectra formation for another pulse and this control 
is effective if the intensity of the control pulse is much larger than that of the investigated 
pulse. While the intensity of the control pulse increases, the squeezing at the frequency 
for which the initial phase was chosen to be optimal, is destroyed. The suppression of the 
quadrature-component fluctuations for the investigated light pulse at higher intensities 
of the control pulse takes place in the frequency domain uj w l/r r . 
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Figure 2. The same as in Fig. |l| but for Qq = 0.5. 



0.5 1.0 1.5 2.0 £2 2 5 3 



0.6 



o 

d 



X 



CO 



0.4 



0.2 









- Shot noise level 


\e 




\d \ 




a 




i.i.i. 


i.i, 



0.5 1.0 1.5 2.0 2.5 3.0 



Figure 3. The same as in Fig. but for fio = 0.7. 



Squeezing spectra of the investigated pulse ( jS5|) for various frequencies (f2 = 0, 0.3, 
0.5) and different intensities of the control pulse as functions of the maximum nonlinear 
phase addition o ,i are shown in Figs. £|, ||, and ||, respectively. One can see from Figs. 
[T|-[| that at higher intensities of the control pulse the spectral density under investigation 
reaches the minimum value in the domain 0o,i — 1 and then the spectral density does 
not practically change. 

The appropriate for our purpose spectra can be also obtained by changing the 
intensity of control pulse and measuring the spectra for defined frequency Q when the 
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Figure 4. Fluctuation spectrum of the quadrature-squeezed component of the 
investigated pulse for the frequency £1 = for the case of optimal initial phase <po,i{t) 
chosen for ilo = and for different intensities of the control pulse 710,2/^0,1 = (a), 
2 (b), 3 (c), 5 (d), 8 (e). Curves are calculated for the time moment t = and 
7i = 72 = 27. 
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Figure 6. The same as in Fig. || but for the frequency = 0.5. 



initial phase of the investigated pulse is chosen to be optimal for various f2 . Such a 
dependence is displayed in Fig. [7[ In some sense, the spectra shown in Figs. f| and [7| 
are similar, and they demonstrate that the measurements carried out at frequencies, for 
which the initial phase was chosen optimal, give almost the same results. From Figs. 
||-[7| one can see that the XPM effect determines the domain of 0o,i values where the 
spectral density of quadrature component fluctuations is independent on the intensity 
of the investigated pulse at the output of nonlinear medium. 



6. Conclusion 



We have developed the quantum theory of two USPs propagation in the nonlinear 
medium with inertial electronic Kerr nonlinearity. In the considered case, the 
propagating pulses are subject simultaneously to the SPM and XPM effects. The 
nonlinear medium is suggested to be lossless and dispersionless, that is, the pulse 
frequencies are off-resonances. Nevertheless, to develop in this case the consistent 
quantum theory of the SPM-XPM effect, it is necessary to take into account a finite 
time of nonlinear response. In the developed approach, we neglect thermal fluctuations 
of nonlinearity, and the causality principle is not violated for the observed values. 

The nonlinear response time defines spectral width of quadrature-squeezed light 
[see Eqs. (|80|), (0)]. The level of suppression of quantum fluctuations depends on the 
nonlinear phase additions due to both the SPM and XPM phenomena. 

It is shown that the frequency, at which the suppressions of fluctuations is 
maximum, can be controlled by adjusting the initial phase of the investigated light 
pulse and the intensity of another pulse. 

The results of this paper can be used in the study of creation of the USPs in 
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Figure 7. Fluctuation spectrum of the quadrature-squeezed component of investigated 
pulse at the frequency — 0.5 for the case of optimal initial phase ipo.i(t) chosen for 
f2o = 0-5 and for different intensities of the control pulse no, 2/^0.1 = (a), 2 (b), 3 
(c), 5 (d), 8 (e). Curves are calculated for the time moment t — and 71 = 72 = 2j. 



nonclassical states and in quantum non-demolition measurements using USPs. 

The approach developed can be applied to generalize the theory of other quantum 
effects associated with the two-mode light propagation in the Kerr media (see references 
in ||) to the case of pulse fields, for example, generation of polarization-squeezed light 
||15|| . This study is in progress. 
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